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We consider the interaction between two rods embedded in a fluctuating surface. The modification 
of fluctuations by the rods leads to an attractive long-range interaction between them. We consider 
fluctuations governed by either surface tension (films) or bending rigidity (membranes). In both 
cases the interaction falls off with the separation of the rods as l/R"^ . The orientational part of 
the interaction is proportional to cos^ [9i + 62] in the former case, and to cos^ [2 {9i + 62)] in the 
^ • latter, where 6\ and 62 are angles between the rods and the line joining them. These interactions 

■ ■ are somewhat reminiscent of dipolar forces and will tend to align collections of such rods into chains. 
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> ■ I. INTRODUCTION AND SUMMARY 

' In addition to their structural role of forming the exterior frames of the cell and its interior organelles and vesicles, 
, lipid bilayers act as the host and regulator of many biophysical and biochemical reactions Inter- and intra- 

cellular recognition and transport, adhesion, regulation of ion concentrations, and energy conversion, are but a few 
of the processes taking place at the membrane. These tasks are carried out by a variety of proteins, glycolipids, and 
other macromolecules that move through the many different lipids that make up the bilayer. The resulting membrane 
is thus far from uniform; there are even examples in which inhomogeneities occur on a larger scale, e.g. domains of 
phase separated lipids or two-dimensional protein assemblies. In modeling the physical properties of the cell, it is 
thus essential to have a good understanding of the interactions between inclusions in fluid membranes. 
. The pursuit of "biologically inspired" materials, which do not posses the full complexity of their natural counterparts, 
(— I ' yet retain some of their useful features, is quite active. Artificial protein assemblies within lipid membranes are now 
Q , routinely produced in the laboratory ||^-|^. Such model- membrane systems have potential applications for targeted 
O ' drug delivery and may also lead to novel applications such as nano-scale pumps, templates, functionalized interfaces, 
and chemical reactors. The appropriate design of such artificial membranes again requires an understanding of how 
inclusions modify the physical properties of the bilayer, and how the membrane in turn contributes to the interactions 
between inclusions. The forces between the inclusions can be broadly subdivided into two categories The first 
^ , category includes interactions that are present in the bulk of the solvent. They include the van der Waals interaction, 
which falls off with separation i? as at long distances. The Coulomb interaction is strongly screened under 

physiological conditions. (Typical ion concentrations are a few hundred millimolar, which give a screening length of 
less than 10 A.) Hydration and structural forces are also short-ranged. The second category includes interactions 
which are mediated by the membrane itself: the inclusion disturbs the lipid bilayer and this disturbance propagates to 
neighboring inclusions (c.f. and references therein). When macroscopic thermal fluctuations are unimportant 

(we refer to this case as T = 0), the resulting interactions tend to be short-ranged. For example, if in the region 
around an inclusion the membrane is forced to deviate from its preferred thickness (~ 40 A), the resulting disturbance 
in the bilayer decays (heals) over a length of order this thickness Two nearby inclusions then feel an interaction 
that falls off exponentially with this characteristic length. 

There are also long-range interactions that are mediated by the membrane. To describe such interactions, it should 
be possible to neglect the microscopic properties of the membrane, and its molecular lipid bilayer structure, and focus 
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on its macroscopic properties. In the long-distance limit, the membrane is well-described by the elastic Hamiltonian 
00, 



n= ds 



-H^ + kK 



(1) 



where dS is the surface area element, and H, K are the mean and Gaussian curvatures respectively. The elastic 
properties of the surface are then described by the tension cr, and the bending rigidities n and R. A finite surface 
tension is in general the strongest coupling in Eq.(|l|) and dominates the bending terms at long wavelengths. This 
term is present for films on a frame, interfaces at short distances, and possibly membranes subject to osmotic pressure 
differences between their interior and exterior. On the other hand, for unstressed vesicles, the surface tension is quite 



small, and may be neglected at wavelengths well below the size of the vesicle [12 -15[]. In this case, the energy cost of 
fluctuations is controlled by the rigidity terms. For simplicity we shall refer to surface tension dominated surfaces as 
films, and to rigidity controlled ones as membranes. 

The long-range interactions between inclusions in a membrane that result from Eq. (|^) were examined in Ref . . If 
the inclusions are asymmetric across the bilayer and impose a local curvature, even at T = 0, there is a long-ranged 
repulsive interaction that falls off with distance as 1/R*. The energy scale of this interaction is set by k and R. On 
the other hand, if thermal fluctuations of the membrane are included (T ^ 0), there is a interaction for generic 

inclusions. The only requirement is that the rigidity of the inclusion differs from that of the ambient membrane 
1^. In particular, the interaction is attractive if the inclusions are stiffer than the membrane. The magnitude of 
this fluctuation-induced interaction is set by fcsT, and is totally independent of the rigidities k and k ( see also 
Appendix ^ . 

In a recent report, we considered the dependence of the fluctuation-induced (T / 0) interaction between rod-like 
inclusions on their orientations Jl6t . The rods are assumed to be sufRciently rigid so that they do not deform coherently 
with the underlying membrane. They can thus only perform rigid translations and rotations while remaining attached 
to the surface. As a result, the fluctuations of the membrane are constrained, having to vanish at the boundaries of 
the rods. Consider the situation depicted in Fig. |l|, with two rods of lengths Li and L2 at a separation R'3> Li. For 
fluctuating films {a ^ 0), there is an attractive fluctuation-induced interaction given by, 

V^{R, eM = cos^ [^1 + 92] + , (2) 

where 9i and 62 are the angles between the rods and the line adjoining their centers, as indicated in Fig. [l[ This 
angular dependence is actually the square of that of a dipole-dipole interaction in two dimensions, with Li and L2 as 
the dipole strengths. The fluctuation-induced interaction on a membrane (cr — 0) is very similar and given by 

Vl,{R, = cos^ [2 (^1 + O2)] + O . (3) 

The orientational dependence is the square of a quadrupole-quadrupole interaction, with the unusual property of being 
minimized for both parallel and perpendicular orientations of the rods. Note that the strength of the interaction is 
the same in both cases. The above fluctuation-induced interactions decay less rapidly at large distances than van der 
Waals forces and may play an important role in aligning asymmetric inclusions in biomembranes. Since orientational 
correlations are often easier to measure than forces, this result may also be useful as a probe of the fluctuation-induced 
interaction. Finally, this interaction could give rise to novel two-dimensional structures for collections of rodlike 
molecules. In particular, the resemblance of the orientational part of the interaction to dipolar forces suggests that a 
suitable way to minimize the energy of a collection of rods is to form them into chains. (If the rods are not colinear, the 
interactions cannot be minimized simultaneously.) Such chain- like structures are observed for ferromagnetic particles 
controlled by similar forces [ pTf . 

In this article we provide the calculational details that lead to the above results. The remainder of the paper is 
organized as follows. Due to its richer complexity, the detailed calculation for the force between rods on membranes 
is presented first in Sec.||. The corresponding calculation for films is described more succinctly in Sec. [II. This is 



followed by a discussion of the results, and comparison with other work in Sec.lV. Various calculational details are 
relegated to the appendices. Asymmetric inclusions have also been recently considered in [^8|, where results similar 
to ours are reported. 



II. MEMBRANES 



We start with a thermally fluctuating planar membrane subject to the Hamiltonian in Eq.(|^). We assume that 
the size of the membrane d, is well below the persistence length ^ m^. In this limit, the membrane undergoes only 
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small fluctuations about a flat state. We may then parametrize the membrane surface with a height function u{r) 
and approximate the full coordinate-invariant Hamiltonian of Eq. (m) by the quadratic form 



(4) 



Since we assume d is large (compared with R and Li), we denote the (finite but large) reference plane by IR^. 

Now consider the situation depicted in Fig. |l|, where two rigid, rod-like objects, are attached to the membrane. We 
shall represent the rods by narrow rectangles of lengths Li and L2, and widths ei and 62; ultimately taking the limit 
of ei — > 0. The rods are constrained to fluctuate with the membrane but, due to their stiffness, can only be tilted or 
translated up and down rigidly. We can parametrize all possible configurations of the rods by 



u{r)\ 



reLi — Oji 



for i = l,2 



(5) 



where we have also used Li to denote the ith rod. To calculate the partition function, we follow a procedure similar to 
Ref. Q and sum over all possible configurations of the membrane, weighted by the corresponding Boltzmann factor, 
and subject to the constraints imposed by the rods via Eq.(||). The constraints may be implemented with the aid of 
delta functions as in Ref. , leading to 
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j Vu{r) Y[ j da^d%, n 6 (u(r') - a, - b, • r') exp 

i—l r'^Li 



Ho 



(6) 



In Eq. (H) we have included only the leading term in an expansion in powers of . As described in Appendix higher 
order terms come from the projection of Li onto the x — y plane, as well as from the integration measure for hi, which 
is on the sphere of unit normals. Since b^ controls the gradient of u(r) at the boundary of Li, the expansion in b.; is 
in the same spirit as the gradient expansion for the Hamiltonian in Eq. (^ . In Appendix ^ we further demonstrate 
that, just as in the case of anharmonic terms that have been neglected in Eq.(^, the higher order terms in b^ left out 
from Eq.(||) are suppressed in the limit d <C ^. Expressing the delta functions as functional integrals over auxiliary 
fields fci(r) defined on the rods, we obtain 



2 

Z = J Vu{r) fl j daA^h j Vh{r) 

d^r{V^u{r)Y + i'Y^ 1 d^nh{ri){u{ri)-ai-h,-Yi) 



(7) 



X exp 



i=l • 
K 



2kBT 

Integrating out u{r), ai, and b^, then gives 



d^nrkiiri) 



Li 



X exp 



keT 
2k 



E 



d'r, / d^rjh{n)G{r^-rj)k,{rj) 

Li JLi 



where 



G(r - r') 



1 

V4 



r - r' I In I r - r' 



(8) 



(9) 



Equation (g) is analogous to the partition function for a pair of plasmas confined to the interior of rods Li and L2- 
The delta functions impose the constraints that the net charge and dipole moments vanish within each rod. When 
the distance R between rods is much bigger than their size (i.e. Li <^ R), we may approximate G(ri — Y2) in Eq.(||) 
by a multipole expansion and keep only the leading term, which comes from the quadrupole moments (see Appendix 



d^r ran ki{r). 



(10) 
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After inserting 



1 = n / dQ«dg(^) exp (q« - / A\ r„ 

i=l-' V ab ^ -^-^i 



n h{r) 



(11) 



into Eq.(^) and performing the multipole expansion, we obtain 

2 = n y"dQ«dgWda,d26, 



X exp 



|-^E dVdV'fc,(r)G(r-r')fc.(/) 



/ dVfci(r) l^fli +bi 



r + r ■ g 



2k ^ 



(12) 



where we have recast the deha functions in Eq.(g) in terms of integrals over a^, and b^. The quadrupole-quadrupole 
interaction 



1 



87ri?2 



-8R • Q(i) • Q(2) . R + 8R • Q(i) • R R • Q(2) . R + 0{1/R^) 



(13) 



(with impHcit summation over repeated a and b) is calculated in Appendix ^ Note that the Green's function in Eq. (|^) 
should also contain homogeneous terms, which reflect the boundary conditions at the outer edge of the membrane, 
r = d. However, we have only used the explicit form of the Green's function in computing the leading terms in the 
multipole expansion. As long as Li and L2 are sufficiently far (compared to R) from the edge, the particular choice 
of boundary conditions at r = d does not modify the leading terms in this expansion. The homogeneous terms can 
therefore be safely suppressed in Eq.(j|). 

We first isolate the integration over fci (r) in Eq. (|l2[) , 



h= y"2?fci(r)cxp|-^ J d\d^r'ki{r)G{r ~r')ki{r') 
d^rfci (r) oi + bi • r + r • g^^^ • r 



(14) 

To perform the above integration, the Green's function in Eq.(^ has to be inverted in the finite region Li. In order 
to do this, we introduce an auxiliary field h(r) and write 



/, = /m(.)exp[-^/^^dV(v^MO)^ 

X JJ^ (5 ih{r') — ai — bi ■ r — r ■ g'"^^ • r 



(15) 



r'eLi 



This is just the partition of a membrane constrained by a single curved rod. After evaluating the contribution on Li 
(via the delta function), we are left with 



/i = exp 



2k 



'Dh{r) exp 



2fcsr JlK^-Li 



d^r {V^hir))' 



(16) 



where the prime indicates that the function h{r)^ and its normal gradient, are constrained to satisfy the following 
conditions on the boundary dLi, of Li, 



Kr)\re^L^ = ai + bi • r + r • g^^^ • r, 



dh{r) 
dn 



9 ( 

— ai + b 
m V 



1 • r + r • g^ 



(17) 
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Now let hoi^r) denote a solution to the biharmonic equation, 

V^/lQ = 

on IR^ — Li with the boundary conditions of Eqs.(p^). We then set 

h{r) =ho{r)+6h{r), 

where both 5h(r), and its normal gradient, are zero on the boundary of Li. Following this change of variables, 



(18) 



Ii — A cxp 



2k r ( (1)'^ 



exp 



2kBT 



(19) 



where 



T>5h{r) exp 



is a normalization constant, independent of a, b, and g, which does not affect the remaining computations. In order 
to solve Eq.(p^) we must specify the boundary conditions at r = d, which are the same as those for u{r). As discussed 
earlier, the results should be independent of this choice, and it is convenient to select 



, I dh 

h\r=d — -TT 

or 



= 0. 



(20) 



As shown in Appendix |^, the solution for the case when the rod is along the y-axis, in the limit d ^ Li, gives 

d\ (V^hoir))^ (21) 



27r ( 



1 



-y [Slbi,+S2biy\ +0(il/d). 



ln(4d/i 

where Si are numerical constants. The second term on the right hand side of Eq.([2l|) is examined in Appendix 
where the irrelevance of higher order terms in b is demonstrated. In the limit d ^ ii, it suffices to keep only the first 
term on the right-hand side of Eq. (|l|) , which gives 



Ii ~ A exp 



kRT 



(22) 



The result of the k2{r) integration in Eq.(p^) is similar, with the index 2 replacing 1, and with the coordinate axis 
appropriately rotated to align with the second rod. The overall expression for the partition function now reads 
(dropping unimportant multiplicative constants) 



2 

Z = Y[ dQWdgWexp 



X exp ■ 



keT 



(23) 



knT 



2k 



where we have set the widths of the rods to zero (i.e. taken the — > limit). The primed indices x' , y' , x" , y" indicate 
that the corresponding components are with respect to the coordinate frames where Li \\ y' and L2 \\ y" ■ We define 
an un-primed coordinate system such that the x-axis is parallel to R and the two rods make angles of 9i and 62 with 
respect to the x-axis as in Fig. |l|. Integration over g yields 

2 

(24) 



knT 



^ = 11 J dQ^'^ S (q« cos^ e, + g« sin 29, + Q,« sin^ 6, 
X S (q« sin^ 9, - Q^l sin 29, + Q« cos^ 9.^ 

Y.^{y^{Qfy-Q^il)sin2e, + Q^^yCos2e^ 



X exp ■ 



2k 
X exp 



keT 
7 

2k. 
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Since we are working in the large- i? limit, the Q integrations are most easily performed by expanding Eq.(p4|) to 
second order in v. After expanding —ksTlogZ, we find the (i?, 0i , 02 )-dependent part of the free energy given in 
Eq.(||). We can re- write this interaction in a coordinate invariant form, in terms of the vector R and the directors Li 
and L2 along the rods as 

2 

+ 0(l/i?'). (25) 



I28 i?4 



2 (2{Li ■ R){L2 ■ R) 



III. FILMS 



We now turn to the analogous calculation for films. Again we use a Gaussian approximation for the Hamiltonian 
in Eq. (^, but keep only the surface tension this time. 



Wo = | / d'r{\7uir)r 



(26) 



All anharmonic corrections to the above expression are unimportant in the limit aa^ ^ ksT, where a is a microscopic 
length. We follow a procedure similar to that described in Sec.(||) but with the differences noted below. The expression 
for the partition function is similar to Eq.(^, with TIq now given by Eq.(p6|). For films however, we cannot justify 
keeping only the leading terms in an expansion in b;. Thus the full rotationally- invariant measure of integration 
on the sphere of slopes hi should be employed (see Appendix Also, the appropriate domain replacing Li is the 

projected length Li = Li/ ^1 + ■ After introducing auxiliary fields ki (r) as in Sec.(||, the analog of Eq.(|) may 

be written as 



X exp 



2 



Vh{r)5{ / d'rMn) 



Li 



2. 



where 



g{v-v') = 



= — In I r — r' 
2tt 



(27) 



(28) 



Note that for films, the dipole moment of ki{r) does not vanish. Expanding Q{v — v') in a multipole expansion and 
keeping only the leading term, which now comes from the dipole moments = Jrfc.i(r), we find 



Z = n/ Vh{r) J dpWrfg.da.— 

' E J- dVd2r'fc,(r)g(r - r')fc.(r') 



X exp ■ 



; J2 j_ d2^fc,(r) [a, + {h, + g,) • r] 



keT 

^ 

2a 



p(l),p(2) 



(29) 



where 



p(l)^p(2) 



7ri?2 



p(l).p(2)-2(p(l).i?)(p(2).7?) 



(30) 



is the effective dipole-dipole interaction, analogous to Eq.(p^), and is the variable conjugate to p*^'^. We integrate 
out ki{r) as in Sec.(||) by introducing an auxiliary field h[r). In this case we must solve a harmonic equation on 
IR^ — Li, instead of the biharmonic Eq.(|l^), with the boundary condition on Li 
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Mr)\redL, = a^ + (b» + g») • r • (31) 

The harmonic problem can be solved either by a method similar to that described in Appendix or by conformal 
mapping. The resulting expression for the partition function reads 



^^^J (1 + bf) 



X exp 



2^3/2 



p(l),p(2) 



(32) 



where the meaning of the primes is the same as in the previous section. One can now see explicitly that the higher 
order terms in the expansion in hi are important in this case. The remaining integrations, except those of biyi and 
62y", can be performed in a straightforward manner. The latter two integrals are rather complicated and in order to 
get a simple result, we restrict to the case cjLF' (^L^/R^) <^ ksT. In this limit, the integrals can be approximated by 
Gaussian forms. After expanding —kBTlogZ, we find Eq.(H), which can be expressed in a coordinate invariant form, 
in terms of the vector R and the directors Li and L2 along the rods, as 



kgT LfLj 
128 R* 



2{Li ■ R)iL2 ■ R) - Li ■ L2 



0(l/i?S). 



(33) 



IV. DISCUSSION 

We shall now discuss some general aspects of the fluctuation-induced interactions in Eqs.(^) and (||). The magnitude 
of the interaction is solely determined by fc^T and is independent of the tension and rigidity coefhcients a and n. 
Thus the effect persists even for rather stiff membranes with k ^ kgT. The only assumption is that the inclusions 
are much more rigid than the embedding surface, thus limiting its fluctuations in their neighborhood. 

For both membranes and films, the interaction falls off with distance as l/R'^. This is a general feature of fluctuation- 
induced forces, including the (finite T) van der Waals interactions, which in d dimensions fall off as 1/R^'^. Since 
the direct van der Waals interactions between inclusions still fall off as 1/R^, the forces mediated through the two- 
dimensional surface will always asymptotically dominate. Of course the dimensional dependence of R'^ is canceled by 
a product of lengths in the numerator. For spherical inclusions, this is given by the product of two inclusion areas 
(see Ref. j|] and Appendix ^) and for rods by the product of the squares of their lengths. Presumably, for general 
shapes, there is a formula that interpolates between these two limits. Another potential extension is to a polymer 
fioating on a membrane. The interplay between the elasticity and shapes of a polymer and membrane, neglecting 
membrane fiuctuations, have been examined in | pO| ; an extension to the case of a fluctuating membrane has also 
recently appeared [ pTf . There is also interesting behavior in the opposite limit of R L for the interaction between 
two parallel semifiexible polymers |2^ . 

Finally, the most interesting aspect of our calculation is the orientational dependence of the force. This is most 
easily discussed for the film, where an intermediate stage involves calculating the angular dependence of a dipole- 
dipole interaction, which is subsequently squared. The final angular dependence is thus that of squared dipolar 
interactions. Similarly, the result for the membrane corresponds to squared quadrupolar interactions. The minimal- 
energy orientations are shown in Fig. ||; note there is a large degeneracy. We also note that these interactions cannot 
be obtained by adding two-body potentials on the rods: To find the orientational dependence of additive forces, let 
us consider an interaction U{\ri — r2\)duidu2, between any two infinitesimal segments of two rods of length L at a 
distance R^ L. Expanding |ri — r2|, and integrating over the two rods, leads to the interaction 

ViR,0i,e2) ^ L'U{R) + ^ (^^^ + U"iR)^ ~ ^ (^^^ ~U"iR)^ (cos20i + cos202). (34) 

Note that the angular dependence is now completely different, and minimized when the two rods are parallel to their 
axis of separation. Presumably both interactions are present for rods of finite thickness; the additive interaction is 
proportional to L^{Le/ R)^ , where e is the thickness. The previously calculated interactions are thus larger by a factor 
proportional to (R/e)^ and should dominate at large separations. 

The unusual dependence on orientation in Eqs.(|2|) and (^ could lead to new types of orientational ordering in 
ensembles of rod-shaped particles. Of course, due to the non-additive nature of the forces, the fluctuation-induced 
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interaction should be calculated separately for each arrangement. However, a cursory examination suggests that 
three-body and higher order interactions fall off with separation as Thus for i? ^ L, a collection of rods can 

be treated as if they interact through additive pair potentials. It is amusing to examine the minimum of such an 
interaction for three rods placed on the vertices of an equilateral triangle. One possible equilibrium configuration is 
a three arm star with the relative angles of 27r/3 between the rods. Interestingly, this so-called "triskelion" structure 
is indeed formed by three rod-like "clathrin" proteins 0. (Another stable configuration has each rod parallel to 
the corresponding side of the equilateral triangle.) Of course, given the relative proximity of the three proteins, it 
is not clearthat the asymptotic interactions of Eq.(||) are applicable to this case. Another generic aspect of dipole 
and quadrupolar interactions is that they are frustrated (i.e. cannot be simultaneously minimized with respect to 
the orientations) if the rod centers are not aligned. There may thus be an overall tendency to arrange rod shaped 
molecules into chains. (Naturally this effect competes with the tendency to aggregate the inclusions together.) We 
hope that the orientational dependent interactions calculated in this paper will provide a fresh perspective on the 
behavior of inclusions in biological membranes. 
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APPENDIX A: THE QUADRUPOLE-QUADRUPOLE INTERACTION 



In this appendix we derive the quadrupole-quadrupole interaction in Eq.([i3|). The starting point is the pairwise 
interaction 



v = 2 dV / d^r' ki(r)G(r -r' - R)k2(r'). 



The Green's function may be written as 



1 



G(r - r' - R) = — I r - r' - R p In I r - r' - R 
Btt 



r2 r'2 



IGtt 



-i?^ 1 + — 7 + — r - 2 



r • r' „ r • R „ r' • R 



i?2 i?2 R2 



- 2- 



R 



R 



(Al) 

(A2) 
(A3) 



In 1 



i?2 R2 R2 



r R r' R 

2 V- 2- 



R 



R 



ln(i?2) 



Since r and r' are limited to the interior of Li and we can expand the above expression in powers of the small 
quantities |r/i?|, |r'/i?|. Because of the constraints 



dVfc,(r) = 0, 



(A4) 



(for i — 1,2) the leading terms in the integration vanish in Eq.(Al). The first non-vanishing term comes from the 
quadrupole moment 



Li 



(A5) 



and is given by Eq.(O). 
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APPENDIX B: INTEGRATION OVER TILT ANGLES 



In this appendix we examine the higher order terms in the tilts of the rods b^, and show that they may be neglected. 
For simplicity we shall focus on bi; similar arguments apply to b2. Whenever possible we drop the index and use 
b = bi and L = Li. The integration for b must be performed over all possible orientations of the rod L in the three 
dimensional embedding space. The manifold of orientations is the unit sphere. In terms of the vector b, defined in 
Eq.(^, the rotation invariant measure on the unit sphere is given by 



dO = 



.2^3/2 



(BI) 



The leading term of the expansion of Eq.(Bl) in b, d?b, was used as the integration measure in Eq.(|^). Additional b 
dependence comes from the projection of the tilted rods onto the x — y reference plane. For example, the conditions 
imposed in Eq.(^) do not really apply for r E L but rather for r in the projected image of L, which is a rod of length 
L/{1 + by)^^'^. Again, in Eq.(||) we have taken the leading order in an expansion in b by setting 1/(1 + by)^/'^ « 1. 

We shall now demonstrate that the higher order terms in b can be neglected (as discussed after Eq.(^). The 
argument is presented explicitly for terms of order b^, but is easily extended to higher orders. Including the first 
corrections to Eq.(|^) results in 



2 

Z = J Vu{r)Y[Vkdr)da,d^b^ (l + T-X+^ybl) 



(B2) 



X exp 



2kBT . 



where and Ty are independent of b. If, as in Sec.|l[ only the leading term is retained, the integration over b leads 
to the constraint that the dipole moment ki(r) must be zero (see Eqs.(^|8|)). Due to the higher order terms in b, this 
constraint is modified, and we have to take into account the dipole moment 



d^r rfci(r). 



(B3) 



Following the same procedure used for the quadrupolc moment in Sec.||, we introduce an auxiliary variable f , via 



1 = y" dpdf exp i{-(^p-J (frrki{r)j . (B4) 
Inserting Eqs.(B4) and (|ll|) into Eq.(p^), and performing the multipole expansion, gives 

Z = I ]JPfc,(r)dQ(*)daidg«db2dbdfdp (1 + T,bl + Tybl) 

X exp ( - ^ / dVd^r'fci (r)G(r - r')fci (/) - ^(b - f ) • p 

-i j d2rfci(r) ai +f -r + r -gfi) -r I X ■•• . (B5) 

In the above equation, • • • denotes factors that are independent of fci(r), b, and f, and identical to the corresponding 
terms in Eq. (p^) with the exception that v [Q'-^-', Q'-^^j is replaced by v [p, Q'-^-', Q^2)j ^ i.e. the multipole energy now 
also depends on p. The integration over fci(r) is the same as in Eq.(|lJ), except that bi is replaced by f. Thus, after 
substituting f for bi in Eq. (pl|) , we are left with the modified integrals 



Z — I dbdfdpexp 



2kBT\n{Ad/L) 



si/J + S2fy) - ^(b - f) • p 



X (l + r,6,2 + Yyby^) {Wo + Wi^pI + Wlypl ^ ■ ' ) 



(B6) 



In the above equation, {Wq, Wix, Wiy, ■ ■ ■} refer to the results of the remaining integrations, which are performed 
after expanding exp[— A:bTz;(p, Q'-^-', Q^2))/2k] in powers of p, and are independent of p, b, and f. After integrating 
over f , and dropping an unimportant constant, we obtain 
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Z = J dbdp cxp [— zb • p] 



(B7) 



X (1 + : 



Tyby 



Wo -Wo 



kBT\TL{Ad/L) pI 



2k 



si 



S2 



+ ■ 



Note that the Wix and Wiy have been dropped since they are subleading in the hmit L. Integrating over b and 
p then gives 



Wo + Wo 



kBT\Ti{Ad/L) fr 



Sl 



tJL 

32 



(B8) 



As discussed in Sec.||, we assume that the size of the membrane is much less than the persistence length ^. Thus, 
the higher order terms in the expansion in Eq.(lB8|) are smaller by powers of 



kBTln{4d/L) _ ln(d/L) 
K ln(^/a) 



<C 1. 



(B9) 



Here we have used the resuh p^ ] ^ « aexp {2'KK/kBT), with a short-distance cutoff a of order molecular size, leading 
to the hierarchy of length scales a < L <^ d To leading order, then, we have Z — Wo, which is independent of 

r^; and Ty and therefore the lowest order term in the expansion in b. It is interesting to note that the above argument 
does not hold for films controlled by surface tension, as discussed in Sec. [II. 



APPENDIX C: SOLUTION OF THE BIHARMONIC EQUATION 

The biharmonic equation (Eq.(|l^)) for a single rod is discussed in detail in this appendix. The problem is to find 
the solution to 

V^h = 0, (CI) 

on a finite disk of radius d from which a line segment of length L along the y-axis has been removed. The boundary 
conditions are 

h{x^Q,-^<y<^)^a + byy + gyyy'', (C2) 



— = 0, -- < y < -) = 5^ + 2g^yy, 



- 2' 
h{d) = 0, 

^h{d) = 0. 
or 

Note that for the boundary conditions, the derivatives are taken before the limit e — > 0. It turns out to be easier to 
impose a weaker boundary condition at r = d, namely 

Md)^o(^), l-m-0[^. (C3) 

Since we have d ^ L, it will suffice to keep the leading terms in the limit L/d ^ 0. Performing the integration by 
parts yields, 

A^r{V^hf (C4) 

IR2-L 



dl V^h- h 



dn &■ 



•n 



d(YR?-L) 

/ dh 
^^Ay [h{0,y)h{y)-—{0,y)f2{y) 



where 
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/i(2/) 



dx 



=0+ 



dx 



(C5) 



f2{y) = v^/.|_„, - v^hl^^_ 



and, as in the text, we have denoted the finite disk of radius d by IR for simpHcity. It is easy to check that the above 
boundary value problem on IR^ — L is completely equivalent to the problem 



y^h = h{y)5{x) + f2{y)^5{x) 



(C6) 



on IR^, provided that the conditions in Eqs.(CS) at r = d are satisfied. The solution to Eq.(C6) can be given in terms 
of the unknown functions fi{y) and f2{y) as 



h{x,y) = 



L/2 

-L/2 
L/2 

-L/2 



dy' Gi{x,y]x' = 0,y')fi{y') 



dy' G2{x, y; x = 0, y')f2{y') 



(C7) 



The Green's functions Gi, which satisfy 

\7^Gi{x,y;x',y') = Six - x') 5{y -y'), 
V''G2{x, y; x', y') = ^S{x - x') 5{y ~ y'), 



(C8) 



dx 



and obey the conditions in Eq.(C3) at r = d, are given by 



Gi{x,y,y') 



G2{x,y,y') 



1 



167r 
1 



x' + iy- y'Y In 



+ {y- y') 



l\21 



1 



an I 2 , 12 \ , ^ ri2 



(C9) 



In 



x'^ + iy- y'f 



, yy 

' d? 



1 - 



rf2 



Note that the boundary conditions in Eq.(|C3|) do not uniquely specify Gi and 6*2, but allow different choices which 
differ by subleading 0{L/d) terms at r = d. Indeed the asymmetry in Gi with respect to interchange of a; and y — y' is 
a result of this freedom. If we require h and dh/dr to vanish at r = d then Gi would be rotationally symmetric. The 
unknown functions /; in the above solution can now be obtained self- consistently by matching to the known forms 
of h and dh/dx on L, as given by the boundary conditions in Eq.(C2). We thus end up with the following integral 
equations 



a + byy + gyyy'^ 



L/2 

-L/2 
L/2 



dy' Gi {x = 0, y; x = 0, y ) fi {y') 



+ 2g,yy = / dy' ^Gs^t = 0, y; x' - 0, y') h{y'). 
J-L/2 dx 



(CIO) 
(Cll) 



(Note that at a; = 0, G2 and dGi/dx are both identically zero.) We start with Eq.(Cll) for f2{y), which is somewhat 
easier to solve. After changing variables to y = (Lcos0)/2 and y' = (Lcos</>')/2, this equation reads 



bx + Lgrcy cos (f> ^ ^ J d(/)'sin0' /2 ^-^cosf/)'^ G'{(f>,(j)') 



where 



G'i 



1 

8^ 



2 ln(2|cos(/i-cos(/>'|) - 2 In 



(C12) 



(C13) 



We now use the expansion |22 
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In (2 1 cos (j) — cos — — ^ — cos ncj) cos n(j)' , 



(C14) 



and define a series 



sin (j)'f2 i ^ cos (j)' I = ^ flm cos mcj)'. 



(C15) 



m=0 



Solving Eq.(C12) for the a^'s gives, to leading order in d L, 

/2 ( ^COS(/i' ) = 



sine 



- 853;^ cos (/)' 



(C16) 



The integral equation for /i requires more care. First, note that the choice of Gi in Eq.(C9) does not lead to a 
vanishing normal derivative at r = c?, unless the condition 



L/2 



d2;V'/i(y') = o, 



-L/2 



is satisfied. Setting up the expansion 



(C17) 



1 cos mo 



m—O 



(C18) 



for /i, this requirement implies 



2bo + &2 = 0. 



The integral equation (CIO) can now be written as 

L'^ L L 

a + Y'^yy + -by cos (j>+ —gyy COS 2<t) = - 



b'smcj^'h [ ^ COS 0')G(0, 



(C19) 



(C20) 



where 



l2 

G{(t>,4>') = — — \{coB(j)- cos (j)'f ln(2|cos(/)-cos(/)'|) 

Ad 



(cos — cos ) In — + 2— - - cos 



1 



COS! 

- COS 2(/) 

oc 

-^COf 

n— 3 



Ad 
L 

2 In 



3 

" 4 

T 



3 1 

-4 + 2 



In 



L2 2 

Ad 

T 



cos 20' 



cos 0' H — cos 3(/)' 
6 



- + In ( ^ ) - i COS 20' + ^ COS 40' 



1 



1 1, /4d 
h - In — 

2 2 V i 

'2 1 



n n — 1 + 1 



1/1 



1 



2 \n+2 n n+1 
1/1 1 2 



2Vn. — 2 n n — 1 



cos(n + 2)0' 
cos(n — 2)0' 



(C21) 



In going to the second form of G{(l),4>') in Eq.( C21 ), we have used the expansion in Eq.(C14) and rearranged the 
resulting expression as a series expansion so that it resembles the LHS of the integral equation. Substituting the 
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expansion of Eq.( Cl^ ) in the integral equation (C20) and equating the coefficients of cosncj) on both sides, we obtain 
the following set of linear equations for the &„, 



{II): 
{III) 
{IV): 



n n — 1 n + 1 



1 / 1 



1 



2 \ n + 2 n n+1 
1/1 1 2 



5 f^d 
- - 2 In — 

2 \ L 



bi + -63 
6 



2 \n — 2 n n— 1 
64, 



1 1 /4d 
-2^2'^[t 



3' 




3 
















4 


.-4^ 







bn+2 

{n>2), 
128 16 



bo - 



1 16 



The solution to the above equations is (to leading order in L), 

bo 
bi 





■4 32 " 




J9yy + 


32 




ln(f)-' 



62 
64 





■4 


32 


m 







= 0, 



384 



(C22) 



(C23) 



and all other 6„ are determined by the recursion relation (I). Putting the results for /i and /2 into Eq.(C4), we find 
Eq.(|l]), with si = S2 = 47r. 



APPENDIX D: SPHERICAL INCLUSIONS 



In this appendix we correct an error in Ref. The expression for Ti. below Eq.(8) in Ref. ||] should read 
{ksT)^ / (1) (1) (1) ^(2)^(2) , r,^(2)A (fcsT)^ 



H 



2k.q 



(Dl) 



This changes the final answer by a factor of 1/2. Thus, Eq.(lO) of Ref. Q for the interaction between two inclusions 
of area A, separated by a distance R, becomes 



V^{R) = -keT 



6A^ 
7r2i?4- 



(D2) 
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FIG. 2. The minimal-energy orientations for two rods in a membrane (a and b) and a film (a only). The energy is minimal 
for all values of 0. 
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